1 = 2 Pion Scattering Length from Two-Pion Wave Function * 
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We present a report on a calculation of scattering length for 1 = 2 S-wave two-pion system from two-pion wave 
function. Calculations are made with an RG-improved action for gluons and improved Wilson action for quarks at 
a' 1 = 1.207(12) GeV on 16 3 x 80, 20 3 x 80 and 24 3 x 80 lattices. We investigate the validity of necessary condition 
for application of Liischer's formula through the wave function. We find that the condition is satisfied for lattice 
volumes L > 3.92 fm for the quark mass range m 2 = 0.273 — 0.736 GeV 2 . We also find that the scattering length 
can be extracted with a smaller statistical error from the wave function than with a time correlation function 
used in previous studies. 



By now the standard procedure employed for 
calculating the scattering length of hadrons from 
lattice QCD is to use the finite-size method pro- 
posed by Luscher, which relates the finite volume 
shift of energy eigenvalues AE to the scattering 
lengths pQ. The formula reads ao/(Ln) = —x — A- 
x 2 -B-x 3 + 0(x 4 ), where x = AE-2m n L 2 /(4Tr) 2 , 
and A = —8.9136 and B — 95.985 are geometri- 
cal constants. The energy shift AE is extracted 
from the asymptotic time behavior of two-pion 
correlation functions. In this way a number of 
calculations has been made for the 1 = 2 S'-wave 
pion scattering length |2I3I4I5I6I7I8I9I1Q|TT] . 

It should be noted that the condition R < L/2 
is assumed for the two-pion interaction range R 
and the lattice volume L 3 in the derivation of 
the Liischer's formula [l|. So far there have been 
studies of the lattice volume dependence of the 
scattering length, but no direct investigation of 
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the interaction range R. It is very important that 
we examine the validity of the necessary condition 
for the Liischer's formula in our current lattice 
simulations. This can be done by investigating 
the two-pion wave functions, which is one of the 
purpose of this article. 

Once one has access to the wave functions, 
one can try to extract the scattering length from 
them. This is the second, and perhaps more in- 
teresting, purpose of this article. Preliminary re- 
sults of the present work was presented at Lat- 
tice'03 H2J. 

Our idea is based on the derivation of the 
Liischer's formula. Luscher found that the two- 
pion wave function 4>(x) on a finite periodic 
box L 3 satisfies an effective Schrodinger equation 
(A + k 2 )(j)(x) = J d 3 y Uk(x, y)<j>{y), where x is the 
relative coordinate of the two pions and k 2 is re- 
lated to the two-pion energy eigenvalue through 
E = 2y/m 2 + k 2 . The function Uk(x,y) is the 
Fourier transform of the modified Bethe-Salpeter 
kernel introduced in Ref. PQ. It is non-local and 
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Figure 1. Wave function on 24 3 lattice at t = 52 
and ml = 0.273 GeV 2 . 

generally depends on energy. 

In the derivation of the formula it is assumed 
that / d 3 y Uk(x,y)<j)(y) ^ only for \x\ < R < 
L/2. In other words the volume has to be suffi- 
ciently large, so that the boundary condition does 
not distort the structure of the two-pion interac- 
tion. Out of this region the wave function satisfies 
the Helmholtz equation (A + k 2 )4>(x) = 0. The 
general solution of Helmholtz equation on a finite 
periodic box L 3 can be written as 



ip-x 



/(P 2 - k 2 ) 



(1) 



up to an overall constant, where T = { p \ p = 
(2ir)/L ■ ft , n £ Z 3 }. Expanding the general so- 
lution JIJ in terms of spherical Bessel ji(x) and 
Neumann ni(x) functions for R < \x\ < L/2, we 
obtain <p(x) — ao(k) ■ jo(k\x\) + Po(k) ■ no(k\x\) + 
• ■ •, where neglected terms are contributions from 
states with angular momentum I > 4. The expan- 
sion coefficients ao(k) and (3o(k) yield the scatter- 
ing phase shift in infinite volume by tan<5n(fc) = 
/3o(fc) / Oio(k) . In particular for the lowest en- 
ergy state of the two-pion, it gives the scattering 
length by (3 (k)/a (k) = a k + 0(k 3 ). The con- 
stant Po(k) is also related to ao(fc) geometrically 
and the relation leads to the Liischer's formula. 

We define the two-pion wave function by 
cj)(x,t) = E R x( 7r ( i? ( :? ) + X,t)n(X,t)S), where 
R is an element of cubic group, and summation 
over R and X projects out the A + sector of the 
cubic group and that of the zero center of mass 
momentum. In order to enhance signals we use a 
source constructed with two wall sources given by 
S = W(t Q )W(t + 1), with W(t) the wall source 
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Figure 2. Ratio U(x) on 24 3 lattice at t = 52. 

at t, in Coulomb gauge fixed configurations. 

We work in quenched lattice QCD employing 
an RG-improved action for gluons at f3 = 2.334 
and an improved Wilson action for quarks at 
Csw = 1-398. The corresponding lattice cutoff is 
estimated as 1/a = 1.207(12) GeV from m p . The 
volumes of lattices (number of configuration) are 
16 3 x80 (1200), 20 3 x80 (1000) and 24 3 x80 (506) 
which correspond to 2.61 3 , 3.26 3 and 3.92 3 fm 3 
in physical units. Quark masses are chosen to be 
ml = 0.273, 0.351, 0.444, 0.588 and 0.736 GeV 2 . 
Quark propagators are solved with the Dirichlct 
boundary condition imposed in the time direction 
and the periodic boundary condition in the space 
directions. The sources are set at to — 12. 

The wave function on a 24 3 lattice at t = 52 
and ml = 0.236 GeV 2 is plotted in Figure □ 
where the horizontal axis is x = \x\. We find that 
the statistical error is very small. 

In order to estimate the two-pion interaction 
range R we construct a ratio U(x) = (A + 
k 2 )(f)(x,t)/(k' 2 (f>(x,t)), where k 2 is obtained from 
the two-pion time correlation function. This ratio 
is expected to vanish out of the two-pion interac- 
tion range R < x. The ratio U (x) on a 24 3 lattice 
at t = 52 for several quark masses are plotted in 
Fig. |21 We find a region U(x) = for x < L/2 for 
all quark masses within the statistical errors. The 
interaction range R tends to be larger for larger 
quark mass. We also find that R ~ 10 (1.6 fm) 
in the worst case of our parameters (the largest 
quark mass, to 2 = 0.736 GeV 2 ). These results 
mean that the necessary condition for the appli- 
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Figure 3. Ratio U(x) on 16 3 lattice at t = 52. 

cation of the Liischer's formula is satisfied on the 
24 3 lattice for all our quark masses. 

The ratio U(x) on a 16 3 lattice at t = 52 for 
several quark masses are shown in Fig. [3] There 
are deviations of U(x) from zero for x < L/2 for 
larger quark masses. These deviation disappear 
for the lightest quark mass, however. This trend 
is also shown for the 20 3 lattice. In the following 
we analyze only data for which we can clearly find 
the region U(x] = for x < L/2. 

We estimate the scattering length by substi- 
tuting the momentum k 2 into the the Liischer's 
formula. The momentum k 2 is obtained by the 
following two methods. In one of them we ex- 
tract the momentum k 2 by a constant fit of the 
ratio V(x) = —A<f)(x,t)/(f>(x,t) for R < x. In the 
other method k 2 is obtained by fitting the wave 
function 4>(x,t) for R < x with the function |JT|I. 

An example of fitting the wave function for a 
24 3 lattice at t = 52 and m 2 . = 0.273 GeV 2 is 
plotted by cross symbols in Fig. ^ The fit works 
very well. 

Finally we compare the scattering length on the 
three lattice volumes in FigQ] We also plot the re- 
sults obtained from the two-pion time correlation 
function on the 16 3 and 20 3 lattices, for which 
we can not clearly find the region U(x) = for 
x < L/2 (data points enclosed by box in Fig. 0J. 
We observe that the scattering length obtained 
from the three methods are consistent within the 
statistical errors. Further the statistical errors 
of those from our new methods are smaller than 
those from the two-pion time correlation function. 
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Figure 4. 1 = 2 pion scattering length ao/m v 
(1/ GeV 2 ) obtained from time correlation function 
(from T), constant fit of V(x) (from V) and fitting 
wave function with Q (from <6). 

We also find no significant volume dependence 
for all quark masses including the data points 
enclosed by the box in Fig. 21 The necessary 
condition of the Liischer's formula is not satisfied 
for these data. However, the effects of deforma- 
tions of the two-pion interaction on the scattering 
length due to finite size effect is apparently small 
compared with the statistical errors. 
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